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Abstract 

The statistics photons in the resonance fluorescence of a qubit excited by microwave and 
radio-frequency (RF) fields have been studied. It has been established that the coherent dissipative 
dynamics of the qubit with allowance for multiphoton emission and absorption of RF photons in 
each act of emission and absorption of a microwave photon results in periodic alternation of photon 
bunching and antibunching. It has been shown that periodically varying statistics photons can be 
smoothly transformed to the purely sub-Poisson statistics by varying the parameters of the RF field. 
The conditions for the formation of the collapse-revival of oscillations in the second-order correlation 
function photons have been determined. The described effects can be implemented on spin and 
superconducting qubits, quantum dots, and qubit-mechanical hybrid systems. 


The measurement of time correlation functions of radiation from light sources in quantum optics 
under lies the determination of its quantum properties and has a long history, including investigation 
of single photon emitters. One of the manifestations of non classical properties of light is photon 
antibunching. This phenomenon was predicted in [1-3] and was observed in resonance fluorescence 
from a monochromatic fieldexcited single atom in a beam [4, 5], a single ion in a trap [6], a single 
molecule in a crystal lattice [7], quantum dots [8], and nitrogenvacancy centers in diamond [9]. The 
characteristic feature of a signal of quantum correlations at quite intense excitation is its modulation 
by the Rabi frequency. Photon antibunching in fluorescence of light scattered by atoms excited by 
intense bichromatic light with nearly resonance frequencies was also studied [10]. 

Quantum properties of an emitter of single microwave photons, including photon antibunching in 
measurements of the secondorder correlation function for intensities, were recently observed with the 
use of artificial single atoms such as superconducting transmons [11, 12]. New features and capabili¬ 
ties of microwave quantum photonics can be manifested in studies of dressed states of qubits at their 
special bichromatic excitation. In this case, the microwave field transfers a qubit from the ground to 
excited state and the resonance frequency of the qubit is modulated by the radiofrequency (RF) field 
with the frequency close to the Rabi frequency of the qubit in the microwave field [13] or much higher 
than this Rabi frequency [14]. At the same time, such studies are of interest because of the possible 
application of such an excitation scheme in quantuminformation technologies [15, 16], to emulate the 
properties of hybrid spin-mechanical systems [17], to measure weak RF fields [18], and in quantum 
amplifiers and attenuators [19]. 

In this work, we study the correlation of photons emitted by spin qubits at their bichromatic 
excitation by the transverse microwave and longitudinal RF fields. 

We consider a twolevel quantum system (e.g., spin qubit) in a microwave field inducing transitions 
between its states and in an RF field modulating its resonance frequency. The Hamiltonian of this 
system can be represented in the form [13, 14] 


H = (1) 

where Hq = is the Hamiltonian of the qubit with the resonance frequency wq; 77j. (t) = 

-I- s~) cosujmwt and 77||(t) = 2uj 2S^ cos{uJrft -|- V”) are the Hamiltonians of the interaction of 
the qubit with the linearly polarized microwave and RF fields, respectively; wi and uj 2 are the cou¬ 
pling constants; tOmw and are the frequencies of the microwave and RF fields, respectively; ip is 


the phase of the RF field; and and are the components of the spin (pseudospin) operator, which 
describe the ground (|g)) and excited (|e)) states of the qubit and satisfy the commutation relations 
[s+,s“] = 2s*, [s*,s^j = ±s^. 

The dynamics of the qubit is described by the master equation for the density matrix p: 

ihdp/dt = [H,p] + iKp (2) 

(below, we set h = 1). Since usually ojilujmw 'C 1, the interaction of the qubit with the microwave 
field is considered in the rotating wave approximation [20]. The superoperator Adescribing decay 
processes is given by the expression [21] Ap = (72i/2)Z)[s“]p + ( 7 i 2 / 2 )U[s+]p + (p/2)U[s*]p, where 
721 and 712 are the rates of transitions from the excited state |e) of the qubit to its ground state ]p) 
and vice versa, respectively; p is the pure dephasing rate; D[0]p = 20p0^ — O^Op — pO^O and O 
is the spin operator. 

We consider two different regimes of the dissipative dynamics of qubits. 

(i) Let Wo — LOmw — rujrf w 0 (r = 0,±1,±2,...) and wi/wr-/ be a small parameter. Then, the 
evolution of the quantum system with time-dependent Hamiltonian o, which is described by mas¬ 
ter equation m, can be considered within nonsecular perturbation theory with the use of Krylov- 
Bogoliubov-Mitropol’skii averaging [14] over fast oscillations after the application of the canonical 
transformation u = exp {—i [wot + ( 2 w 2 /Wf./) sin(wr/t -I- '*/’)]} fo Eq. ([2]). After the implementation of 
both procedures, the equation for the modified density matrix which describes the rth resonance, 
can be obtained from master equation Q in the form idprjdt = [H^f f(r), Pr] + iAp^.. Here, the 
time-independent effective Hamiltonian Heff{r) describes r-quantum transitions. As a result, 

H ^ Heffir) = Ho{r) + ( 3 ) 


where 


— (^0 ^mw TUJrf')S , 


= iH(r) (s+ -b s ) , 
= ABs{r)s^, 


ft{r) = UJiJ-riz), 


Assir) = 


2 _ n)ujrf 

n^r ■' 


Jniz), 


H(r) = uiiJ-r{z) is the rf-field-renormalized Rabi frequency of the qubit in the microwave field, 
ABs{r) is the Bloch-Siegert shift of the frequency of the resonance transition of the qubit (which is 
nonzero only for r 7 ^ 0), and 2 = 2 w 2 /wr/ is the argument of the Bessel functions (which determines 
the magnitude of the effective interaction il{r) with the bichromatic field). Taking into account the 
Bloch-Siegert shift, we redefine the resonance conditions as wq + ABsif) ~ <^mw — rtOrf ~ 0. The 
solution of the master equation for the density matrix of the qubit p^ in the reference frame rotating 
at the frequency ujmw can be obtained in the explicit form (it is assumed that the qubit at the initial 
time t = 0 was in the ground state ]p)): 


‘W = - 


721 — 712 


X { 2(7 + a)s"' + in{r) [s~'~ fr{t) - H.c.] } + 


(4) 


R(7,i-7i2) 

V 4Ha(r) _ flair)-ij 


X [{H(r) [s+fr{t) 


H.C.] — 2 [nc((r') + ia] s^} + iL.c.] 


where frit) = exp{-i[rip + zsiniujrft + if)]}, flair) = \/fl'^ir) - o?, a = ( 7 _l - 7 ||)/ 2 , 7 = 
( 7 ± + 7||)/2 and 7 || = 712 -b 721 and 7 _l = ( 7 || -b p )/2 are the rates of energy and phase relaxations of 
the qubit, respectively (below, we set 712 ~ 0 ). 

From Eq. o, we obtain the following expression for the population of the excited level of the 
qubit absorbing one microwave photon with simultaneous absorption or emission of r RF photons: 

1 H^(r) 

2 9?ir) -b7||7_L 
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With the use of Eq. O and quantum regression theorem [20] , an expression can be obtained for cor¬ 
relation functions of field amplitudes proportional to the spin operators of the qubit. The normalized 
second-order correlation function for the rth resonance can be represented in the form 


9 ?Ht) = (e|p7‘(t)|e) |p(o)=| 3 )( 3 | /{e\p7\t oo)|e), 


( 6 ) 


where {e\p!^°*{t)\e) |p(o)=|g)(g| means that the solution of the master equation for the density matrix 
should be sought for the situation where the qubit at the initial time is in the ground state \g). 
According to Eqs. (O and ([B]), oscillations of the correlation function are caused by coherent 

oscillations of the population of excited levels of the emitter. For r = 0, correlation is absent; i.e., 
gr TO) = 0. This indicates the quantum nature of radiation leading to the photon antibunching effect. 
In the case under consideration of the excitation of the qubit by the bichromatic field, emission and 
reemission are complicated. The modulating RE field modifies the energy spectrum of the system: 
the ground and excited states of the qubit are supplemented by levels splittings between which are 
specified by the frequency ujrf (Fig. la) [14]. Owing to the modified spectrum of the qubit at the 
main resonance (wq ~ >-^mw = 0), multiphoton absorption of a certain number of RE photons and 



Fig. 1. (a) Scheme of multiphoton transitions excited by bichromatic radiation in the case r = 0. 
Thick lines correspond to the levels of the initial spin system that is not modihed by radiation, (b) 
gr (t) for r = 0,—1 at bichromatic excitation of qubit, comw = ojq, = 0.1 MHz, a;r-//27r = 1 

MHz, z = 0.5 z = 2.0, 7 J. = 0.05 MHz, 7 || = 0.01 MHz. 




























emission of the same number of photons occur simultaneously with the emission or absorption of a 
microwave photon. These processes are taken into account in the zeroth order Bessel function and 
reduce the probability of absorption of a microwave photon owing to the renormalization of the Rabi 
frequency: wi —>■ 17(0) = uJiJo(z), z = ‘IfjJijiOrf- With an increase in the parameter z, the function 
Jo(z) decreases and is zero at z\ k, 2.41, zi k, 5.52, etc. Consequently, with an increase in z, the effec¬ 
tive frequency of oscillations ila{0) = -v/wj Jq (x) — decreases (if a;iJo( 0 ) > a), becomes zero at a 
certain z value when Jo(z) = a, and then becomes imaginary and the oscillation-relaxation behavior 
is replaced by the relaxation behavior (see Figs, lb and 2a). Therefore, purely sub-Poisson statistics 
are established instead of periodically varying statistics of photons in a time window comparable to 
the relaxation time I/ 7 . 
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Fig. 2. (Color online) Correlation function gr Rr) at bichromatic excitation of the qubit 
versusz = 2 w 2 /wr/ for ojq + ABs(r) — ujmw — roJrf ~ 0, uJil2Tr = 0.1 MHz, 7 j_ = 0.05 MHz, 7 || = 0.01 
MHz. (a) r = 0. (b) r = —1. The dotted (z = 0.5) and dashed (z = 2.0) lines indicate the positions 
of the signals shown in Fig. lb. 


We note that the role of the function Jo{z), which appears owing to multiphoton absorption and 
emission of the same number of RF photons, is similar to the role of the Debye-Waller factor in solid 
state physics. This factor, e.g., takes into account the elastic scattering of phonons at the interaction 
of the crystal lattice with electrons, which leads to the formation of phononless lines in optical spectra. 
In this respect, oscillations of the nanomechanical system interacting with the qubit can serve as the 
RF field; as a result, the properties of, e.g., the hybrid spin-mechanical system can be emulated by 
means of the modulating RF field. 

At side resonances, wq -I- AB5(r) — uJmw T 1^1 ujrf = 0, the absorption of a microwave photon is accom¬ 
panied by the absorption (-) or emission (-I-) of RF photons whose number is larger than the number 
of emitted or absorbed RF photons by |r|. In this case, the frequency of oscillations of the function 
gi^\T) is given by the expression flair) = ^/f^^ir) — a^. The implementation of both the main and 
side resonances is accompanied by multiphoton absorption and emission of the same number of RF 
photons, which is taken into account in the Bessel functions in the expressions for renormalized Rabi 
frequencies, D(r) = ujiJ-riz). 

Variation of these frequencies owing to the variation of the parameters of the RF field (z value) is 
manifested in the oscillation behavior of the correlation function g). Rr), which vanishes in the limit 
r —>■ 0 (Figs, lb and 2b). As is seen, in contrast to the main resonance, an increase in z for the side 
resonance results in the transition from the sub-Poisson statistics of photons to periodically varying 








statistics. Thus, in contrast to the monochromatic excitation, photons are emitted and absorbed 
cascade-by-cascade rather than one-by-one. As a result, the experimentally observed emission and 
further absorption of only one microwave photon at the main resonance and for the case of side res¬ 
onances, e.g., for r = — 1, the emission and absorption of one microwave photon and one RF photon 
occur with smaller probabilities because of the corresponding factors Jq{z) and , taking into 

account directly unobserved processes of multiphoton absorption and emission of RF photons. 

(ii) We assume that >> W 2 , a; 2 /wr/is a small parameter, the microwave field is strong 
(uif ^ 7||7_l), and the conditions of the exact resonance of the microwave field with the qubit 
(ojjnw = Wo) and Rabi resonance (wi = uj^f) are satisfied. Then, the density matrix of the qubit 
in the bichromatic field can be represented in the rotating reference frame in the form 

Prot{t) = 2 “ 4 (7) 




— s +2s^)—H.c. -|--sin^ 


X (R{t) { 

+R * (t) { 


^ 2s") - 7J.C.] + 2(s+ + s-)} 


_ g- + 2s") - H.c] - 2(s+ -b . 


where R[t) = e * [cosVLpt — sin , flp = \/^\ — / 3 ^ , 7 = (7_l + 7||)/2, l' = 1 — P and 
P = 7||/4. Using the regression theorem [20] and Eq. Q, we obtain the normalized second-order 


correlation function 


5^^^(r) = l —e cosi/;cos(a;ir-b ^/’) 


( 8 ) 


— - sin^e ^ * (sin [(wi -b -b V'j + sin [(wi — 0 ^)t -b ip\ 

+ 7 ^ [cos [(wi -b n/ 3 )r -b V”] - cos [(wi - Vtp)T + p]] 

According to Eq. dH), g^\o) = 0; i.e., the radiation of the qubit is of a quantum character. The mi¬ 
crowave field at the bichromatic excitation induces Rabi oscillations with the frequency wi and splits 
each of the levels of the qubit into two sublevels with the energy gap , whereas the RF field induces 
Rabi oscillations with the frequency u )2 and splits each of the formed sublevels into two sublevels with 
the energy gap UJ 2 - Figure 3a illustrates the double splitting of the upper level of the qubit in the 
bichromatic field. As a result, the emission spectrum of the (qubit + bichromatic field) system includes 
three Mollow triplets: uJmw, Wmu)±W 2 ; Wmuj+wi, u)raw+oJikoJ 2 andwmu; —wi, uJmw— 0 Jiku ]2 [16]. Os¬ 
cillatory terms in Eq. ([8|) for the correlation function (r) correspond to quantum transitions with 
the frequencies wi and wi ±a; 2 . As is shown in Eig. 3b, the collapse-revival effect is clearly manifested 
in oscillations of g^^\T) ; i.e., the intensity of radiation of the qubit vanishes periodically at the times 
Tn = (2n + l)7r/n,g and then is restored. This occurs because oscillations of the population difference 
between the qubit levels at the frequency wi are modulated by rf-field-induced Rabi oscillations with 
the frequency Change in the amplitude of the RE field can retard or hasten the development of 
collapse and revival of oscillations in the correlation function, thus changing the statistics of emitted 
photons (Fig. 4). The oscillatory behavior of g^'^'^ir), including the collapse-revival effect, is sensitive 
to the phase of the RF field (see Eq. (8) and Fig. 3b). At ip = 0, the effect of the RF field in-phase 
with the microwave field is completely eliminated and oscillations occur only at the frequency wi. If 
the phase of the RF field is shifted by 7r/2 with respect to the phase of the microwave field, oscillations 
at the frequency wi disappear and oscillations appear at the sum (wi +ilp) and difference (wi —^p) 
frequencies, leading to the collapse-revival effect. At a random phase of the RF field, Eq. (|S|) is 
averaged over the uniform distribution of ip from 0 to 27r. In this case, oscillations at three frequencies 
(wi, wi ± ilp) make comparable contributions and the collapse-revival effect is most pronounced. In 
the laboratory reference frame, groups of microwave and RE photons at the frequencies corresponding 
to three Mollow triplets are involved in emission and absorption. If the RF field is absent, Eq. (|8|) is 
reduced to the known expression [3, 20] for the case of monochromatic excitation when the correlation 
function oscillates with the frequency wi. 

Thus, the second-order correlation function for resonance fluorescence of the qubit excited by the 
microwave and RF fields has the properties confirming the emitted photon antibunching effect. It was 
established that emission and absorption acts are multiphoton. At wi <gC Wj./, multiphoton emission 
and absorption of RF photons occur simultaneously with radiation or absorption of each microwave 










photon. These processes can be controlled by the RF field and by varying the behavior of the correla¬ 
tion function until the complete suppression of oscillations and transition to the antibunching regime 
determined by the dephasing rate. At wi = cOrf, there is the effect of collapse-revival of oscillations 
of the correlation function caused by R-field-induced Rabi oscillations between dressed states of the 
qubit. The predicted properties of the correlation function can be potentially applied not only for 
spin qubits but also for artificial single atoms, such as superconducting qubits or quantum dots, as 
well as for qubit-mechanical hybrid systems. 



Fig. 3. (a) Scheme of transitions excited by bichromatic radiation in the case ojrf ~ wi in the 
rotating reference frame, (b) Correlation function at uimw = wq, uJrf = wi, uji/2tt = 1 MHz, 

u}2/2tt= 0.1 MHz, 7_L = 0.05 MHz, 7|| = 0.01 MHz. 





MHz, 7 j_ = 0.05 MHz, 7 || = 0.01 MHz, ijj is a random variable. The dashed line corresponds to the 
signal presented in Fig. 3b. 
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